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Abstract 

The local Gauss law of quantum chromodynamics on a finite lattice is investi- 
gated. It is shown that it implies a gauge invariant, additive law giving rise to a 
gauge invariant ^-valued global charge in QCD. The total charge contained in a 
region of the lattice is equal to the flux through its boundary of a certain ^-valued, 
additive quantity. Implications for continuous QCD are discussed. 
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1 Introduction 



Quantum chromodynamics (QCD) is one of the basic building blocks of the standard 
model for describing elementary particle interactions. During the last decades QCD was 
quite successful, e.g. in describing deep inelastic scattering processes within the framework 
of perturbation theory at one hand and "measuring" certain types of observables using 
(nonperturbative) lattice approximation techniques on the other hand. Nonetheless, we 
are still facing the basic challenge, which consists in constructing an effective microscopic 
theory of interacting hadrons out of this gauge theory. For this purpose, nonperturbative 
methods for describing the low energy regime should be developed. In particular, the 
observable algebra and the superselection structure of this theory should be investigated. 
The present paper is a step in this direction. 

We stress that standard methods from algebraic quantum field theory for models which 
do not contain massless particles, see 0, do not apply here. Some progress towards an 
implementation of similar ideas for theories with massless particles has been made, for 
the case of quantum electrodynamics (QED) see || and M and further references therein. 
In particular, Buchholz developed the concept of so called charge classes and found a 
criterion for distinguishing between the electric charge (carried by massive particles) and 
additional superselection sectors corresponding to inequivalent asymptotic infrared clouds 
of photons. For some attempts to deal with the nonabelian case we refer to papers by 
Strocchi and Wightman, see || and ||. 

In QED, the notion of global (electric) charge is easy to understand. This is due 
to the fact, that in this theory we have a local Gauss law, which is built from gauge 
invariant operators and which is linear. Thus, one can "sum up" the local Gauss laws 
over all points of a given (spacelike) hyperplane in space time yielding the following gauge 
invariant conservation law: The global electric charge is equal to the electric flux through 
a 2-sphere at infinity. On the contrary, in QCD the local Gauss law is neither built from 
gauge invariant operators nor is it linear. The main point of the present paper is to show 
that it is possible to extract from the local Gauss equation of QCD a gauge invariant, 
additive law for operators with eigenvalues in Z 3 , the center of 577(3). This implies - as 
in QED - a gauge invariant conservation law: The global Z 3 -valued charge is equal to a 
Z 3 -valued gauge invariant quantity obtained from the color electric flux at infinity. There 
is a paper by Fredenhagen and Marcu, see f7|, dealing with a Z 2 -gauge theory (with Z 2 - 
valued matter fields) on the lattice. These authors were able to construct the ground state 
and charged states of this model. For some regions in the space of coupling constants, the 
thermodynamic limit of charged states was controlled. It is likely that similar techniques 
will be relevant for studying Z 3 -valued colour charged states of QCD. 

We stress that our main observation is universal in the following sense: Suppose some- 
body had constructed a nonperturbative version of QCD rigorously. Then our construc- 
tion of the global charge would apply. Since, however, such a theory is not at our disposal, 
for a rigorous discussion we have to restrict ourselves to the case of QCD approximated 
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on a finite lattice. For basic notions concerning lattice gauge theories (including fermions) 
we refer to |IJ and references therein. Thus, we consider QCD on a finite, regular cubic 
lattice in the Hamiltonian approach. Our starting point is the notion of the algebra of 
field operators. Since, up to our knowledge, this has never been published, we discuss 
it in some detail. By imposing the local Gauss law and gauge invariance we obtain the 
observable algebra 0(A) for QCD on the lattice. For an analogous investigation of the 
superselection structure of QED on the lattice we refer to || and for a (rather heuristic) 
discussion of continuum QCD within the functional integral approach (using only gauge 
invariant quantities) we refer to @. For QED, the mathematical structure of O(A) has 
been completely clarified, see [|10j. A similar investigation for QCD seems to be a very 
difficult problem, which will be addressed in the future. In this paper, we concentrate 
on a careful analysis of the local Gauss law and its consequences for the notion of global 
charge in QCD. 

Our paper is organized as follows: In Sections |2| and ^ we discuss the algebra of fields 
and the observable algebra for QCD on a finite lattice. In Section ^ we analyse the local 
Gauss law and invent the gauge invariant notion of the local charge densitiy. Using this 
notion, in Section [5] the global charge is defined and the flux law for QCD is discussed. 
Here also some heuristic remarks concerning the continuum case are added. 

2 The Field Algebra for Lattice QCD 

We consider QCD in the Hamiltonian framework on a finite, regular 3-dimensional lattice 
A. We denote the set of oriented, i-dimensional elements of A by A*, i = 0,1,2,3. Such 
elements are (in increasing order of i) called sites, links, plaquettes and cubes. The set 
of non-oriented z-dimensional elements will be denoted by |A|\ If, for instance, (x,y) G 
A 1 is an oriented link, then by \(x, y)\ G | vV | 1 we mean the corresponding non-oriented 
link. Instead of using a concrete Hilbert space representation (e.g. the Schrodinger 
representation), we give an abstract definition of the field algebra in terms of generators 
and defining relations. 

The basic fields of lattice QCD are quarks living at lattice sites and gluons living on 
links. The field algebra is thus, by definition, the tensor product of local fermionic and 
bosonic algebras: 

^(A):= (g) T x (g) T\ {x , y)] . (2.1) 

xg|A|Q Kx^lelAl 1 

We impose locality of the lattice quantum fields by postulating that the local algebras 
corresponding to different elements of A commute with each other. 

The fermionic field algebra T x associated with a lattice site x is the enveloping algebra 
generated from the complex Lie super-algebra C x of canonical anticommutation relations 
of quarks. In terms of coordinates, the quark field is given by 

\A\°3x^(iP aA (x))eC x , (2.2) 
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where a stands for bispinorial and (possibly) flavour degrees of freedom and A = 1,2,3 
is the colour index corresponding to the fundamental representation of SU(3). The con- 
jugate quark field is denoted by ip* a A (x) , where we raise and lower indices by the help 
of the canonical hermitian metric tensor g A B in C 3 . Finally, the (nontrivial) canonical 
anticommutation relations for elements of C x read: 

[r a A (x)^ bB (x)] + = 5 B A 5 ab . (2.3) 

Passing to self-adjoint generators (real and imaginary parts of ip), we observe that T x is 
a Clifford algebra. 

The bosonic field algebra J~\( x ,y)\ associated with the non-oriented link \{x, y)\ will be 
constructed in terms of its equivalent copies J-( x , y ) and J-( y>x ), corresponding to the two 
orientations of the link (x, y) . We will see that there is a natural identification of these two 
algebras, induced from the vector character of both the gluonic potential A and the colour 
electric field E of the underlying continuum theory: under the change of orientation of 
the fc'th axis both Ak and E h change their signs. 

The bosonic field algebra Ffay) is, by definition, the tensor product 

F(x,y) '■= £(x,y) <8> A( x ,y) , (2-4) 

where £{ x , y ) is the enveloping algebra of the (real) Lie algebra S( x ,y) — su(3) and A( x>y ) = 
(7°° (577(3)) is the commutative *-algebra of smooth functions on the Lie group SU(3) . We 
note that the above tensor product is naturally endowed with the structure of a crossed 
product of Hopf algebras, given by the action of generators of £( x>y ) on functions, see 
flTTf . We identify the tensor product of elements of £( x , y ) and functions with their product 
as differential operators on SU (3) . This way, Ffay) gets identified with the algebra of 
differential operators on the group manifold. For our purposes, we endow T( x ^ y ) with 
a Lie algebra structure. Thus, we have to define the commutator between a generator 
e G £(x,y) of £{x, y ) and a function / G A/ X>y ) , 

[eJ]:=e R (f), (2.5) 

where e (/) denotes the derivative of / with respect to the right-invariant vector field e R 
generated by e , 

e R (f)(g) ■= -f- /(exp(se) • g) , g E SU(3) . (2.6) 
as s=o 

Now we give an explicit description of J-{x,y) in terms of generators and defining rela- 
tions. The algebra A( X} y) is generated by matrix elements of the gluonic gauge potential 
on the link (x, y) , 

A 1 3 (x, y) -> U A B (x, y) G A [x , y) , (2.7) 

where A, B = 1, 2, 3 are colour indices. Being functions on 577(3) , these generators have 
to fulfil the following relations: 

(U A B (x,y))*U A c (x,y) = 6 C B 1 , (2.8) 
e A BcU A D (x,y)U B E (x,y)U c F (x,y) = e DEF 1 . (2.9) 
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The algebra £( x , y ) is generated by colour electric fields, spanning the Lie algebra £( x ,y)- 
Choosing a basis , i = 1, . . . , 8, of su(3) we denote by {Ei(x,y)} the corresponding 
basis of £( x ,y), 

A 1 3 (x,y) -> :=t« G £( X)1/ ) . (2.10) 

These generators are self-adjoint (real), i£* = . In the sequel we take as the basis the 
hermitean, traceless Gell-Mann matrices t A B , normalized as follows: 



Y, t * A Bti a D = 6£S$-±6iS°, (2.11) 



or, equivalently, £ f B t 3 A = 5^. We will also use the following traceless matrix built 
from these fields 

E A B (x,y) = J2^y)U A B ■ (2.12) 

i 

Since the coefficients ti A B are complex, these fields are no longer self-adjoint: 

(E A B (x,y))* = E B A (x,y). (2.13) 

The commutation relations between elements of su(3) translated to the language of these 
fields read 

[E A B (x, y), E c D (x, y)\ = 5 c B E A D {x, y) - 5 A D E c B (x, y) , (2.14) 

whereas the commutation relations (|2.5|) between Lie algebra elements and functions, 
rewritten in terms of generators, take the following form: 

[E A B (x, y), U c D (x, y)} = 5 c B U A D {x, y) - ±5 A B U c D (x, y) . (2.15) 

Observe that, for every link (x,y) , we have a model of quantum mechanics with 
configuration space being the group manifold SU(3) . The matrix elements of the gluonic 
potential play the role of functions of position variables, whereas the colour electric fields 
play the role of (non-commuting) canonically conjugate momenta. Formula (|2.15 ) is the 
analog of the canonical commutation relation [p, q] — — i . 

Now, we describe the transformation law of these objects under the change of the link 
orientation. This law is an isomorphism between bosonic field algebras: 

Z(x,y) ■ F(x,y) -» F(y,x) ■ (2.16) 

The vector character of the continuum gauge potential implies that the (classical) SU(3)- 
valued parallel transporter g(x,y) on (x,y) transforms to 5 f ~ 1 (x, y) under the change of 
orientation. Thus, the change of orientation on (x, y) induces the transformation 

SU{3) 3 g - i(g) = g' 1 6 SU(3) (2.17) 
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on configuration space. This transformation lifts naturally to the field algebra F(x,y) ■ 
Indeed, for elements of Af x>y ) it implies the following transformation law 

•A(x, y ) 3 f — > f G A(y, x ) , (2.18) 

where f(g) := fig' 1 ) ■ We thus put I( Xjy )(Q, f) := (0, /). It may be easily proved that the 
unique *-isomorphism, which fulfils this requirement is given by 

Z {x>y) (e,f):=(e L J) , (2.19) 

where e L is the left invariant vector field on 577(3) generated by — e. It acts on functions 
in the following way: 

e L (J)(!l) = T f(9 ■ exp(- S e)) . (2.20) 
as s=o 

Observe that e L is not an element of £( y , x ), but it can be expanded with respect to 
right-invariant vector fields with coefficients being functions on the group. Actually, the 
following formula is easily proved: 

8 

eL = -X>(<? e <rW , (2.2i) 

where £, is any orthonormal basis of su(3). 

To prove that It x ,y) is an isomorphism indeed, we observe that it is obviously a bijective 
mapping between generators. It extends to an isomorphism of field algebras if we prove 
that it preserves the commutator. This is obvious for elements, which either belong both 
to A( X)V ) or to £( x ,y)- Thus, it is sufficient to consider the commutator of a Lie algebra 
element e with a function / . Using the identity 



d 



ds 



d 



s=o ds 



e R (f)(g) = -r /(exp(ae) • g) = ~r fig' 1 ■ exp(-ae)) = e L {f){g- 1 ) , (2.22) 



s=0 



and applying I( x , y ) to the functions on both sides, we obtain due to ( |2.5| ) and (|2.18 



T {XtV) ([e,/]) = l [x , y) {e R {f)) = e\f) = [I ix>y) (e),l {x>y) (f)] . (2.23) 

This shows that 1( x , y ) is an isomorphism, indeed. We add that, in more abstract terms, 
it is given by I( x , y )(e, f) = (-i*(e),i*(/)) . 

Field configurations, which are related under I( x , y ) will be identified as different rep- 
resentations of the same object. Thus, the bosonic field algebra J~\{ x , y )\ , associated with 
the non-oriented link | (x, y) \ , is defined as the subalgebra of J-( x ,y) x J~(y,x) , obtained by 
this identification: 

^\(x,y)\ ■= {{l(x,v)ihv,x)) 6 F{x,y) X F{y,x) =Z{x,y){l{x,y)) } ■ (2.24) 
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Projection onto the first (resp. second) component gives us an isomorphism of J-u X; y)\ with 

F(x,y) ( res P- F(y,x))- 

The above transformation law yields relations between generators of the two algebras. 



For functions on SU(3), formula ( |2.8| ) enables us to rewrite (|2.18|) as follows: 

U A B (y,x) = (U B A (x,y)y . 
For Lie algebra elements, formula (]2.21|) applied to e = U reads: 



(2.25) 



Efax) = -J2 uA B(x,y)tj B C U C ' D (y,x)ti D 'a Ei(x,y) 



(2.26) 



or, in terms of generators E A B (x,y) , 

E A B (y,x) = -U A D (y,x)U c B (x,y)E D c (x,y) . 

To summarize, the field algebra J 7 (A) is the *-algebra, generated by elements 

{r A (x),U A B (x,y),E A B (x,y)} , 



(2.27) 



fulfilling relations (Q, (pj|) , ( gTTg) , (ggp and (|2~27| ), together with canonical (anti-) 
commutation relations 



ir a A (x)^ bB (y)] + 
[E A B (x,y),E c D (u,z)} 

[E A B (x,y),U c D (u,z)} 



5(x-y)5 B A 5 ab , 

5(x - u)8(y - z) (5 c B E A D (x, y) - 8 A D E c B (x, y)) , 
+ 6{x - u)S{y - z) U c B U A D (x, y) - ^8 A B U c D (x, y) 

- 5(x - z)5{y - u) (5 A D U c B (y,x) - U A B U c D (y,x) 



3 The Observable Algebra 



The observable algebra O(A) is defined by imposing the local Gauss law and gauge in- 
variance. 

The group of local gauge transformations acts on ^(A) by automorphisms x — > g A B (x) 
as follows: 



r A (x) -> g A B (x)r B (x), 

U A B (x,y) -> ^cix^Dix^ia^fM 



E A B (x,y) -> / c (x)E c D (x,i/)((7- 1 ) D B (x) . 



(3.28) 
(3.29) 
(3.30) 
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It is easy to check that these transformations are generated by 

Q A B {x) := p A B {x) - £ E A b(x, y) , (3.31) 

y 

where 

p a b (x) = E U* aA (x)r B (x) - y A r aC (x)rc(x)) (3.32) 

is the local matter charge density. Observe that p A A (x) = . 

To implement gauge invariance we have to take those elements of J- (A) , which com- 
mute with all generators Q A B {x) ■ Thus, the subalgebra of gauge invariant fields is, by 
definition, the commutant A(Q)' of the algebra A{Q) , generated by the set {G A b(%)} ■ 

The local Gauss law at x G A has the form 

Y J E A B (x,y) = p A B (x) , (3.33) 
v 

with the sum taken over all points y adjacent to x . Imposing it on the subalgebra A{Q)' 
of gauge invariant fields means factorizing the latter with respect to X(A) D A(G)' , where 
X(A) is the ideal generated by ( |3.33| ). Thus, the observable algebra is defined as follows: 

0(A) : = A(gy/{i(A) n A(gy} . (3.34) 

Obviously, every element of 0(A) is represented by a gauge invariant element of ^(A) 
and 0(A) can be viewed as a *-subalgebra of .F(A) generated by gauge invariant bosonic 
combinations of U and E and by gauge invariant combinations of ip and ip* of mesonic 
and baryonic type, with the Gauss law inducing some identities between those generators. 

It will be shown in the next section that there is a (gauge invariant), additive law, 
obtained by combining these two equations, which characterizes the local colour charge 
density carried by the lattice quantum fields. Additivity will allow us to obtain the global 
colour charge by adding up the local Gauss laws. In a subsequent paper we will show 
that the irreducible representations of 0(A) are labeled by this global charge. This way, 
we get the physical Hilbert space 'H phys as a direct sum of colour charge superselection 
sectors. The similar problem for QED on the lattice was solved in ||. 



4 The Local Charge density 



In this section we are going to analyze the local Gauss law (|3.33|) . Suppose, for that 
purpose, that we are given a collection of operators F A B in a Hilbert space 7i, fulfilling 
F A a = and (F A B )* = F B A } realizing the canonical commutation relations for the Lie 
algebra su(3): 



[F' 



B ■ 



D 



S C nF A n -8 A nF 



C 



B 



(4.1) 
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The field algebra .F(A) provides us with two basic examples of this type: the electric 
field E A B (x,y) on each lattice link, see (|2.14|) , and the charge operator p A B (x) at each 



lattice site, given by formula (|3.32| ). Indeed, due to canonical anticommutation relations 
( |2.3| ), the latter fulfills (|4.1|) . Thus, the operators occurring on both sides of the local 
Gauss law fulfill (|4.1| ). 



Throughout this paper, we assume integrability of the Lie algebra representations 
under consideration. This means that for each F there exists a unitary representation of 
the group SU (3) 

SU(3) 3 g - F(g) E B{H) , (4.2) 

associated with F. 

It is easy to check that if F and G are two commuting representations of su(3) then 
also F + G is. Indeed, if F(g) and G(g) are representations of SU(3) correspond- 
ing to F and G, then F + G may be obtained by differentiating the representation 
SU(3) 3 g — > F(g)G(g) E B(H), where B(TC) denotes the C*-algebra of bounded opera- 
tors on 7i. Moreover, — F* is also a representation of su(3), corresponding to the following 
representation of SU(3): SU(3) 3 g^ (Fig- 1 ))* E B(H). 

Such a collection of operators is an operator domain in the sense of Woronowicz (see 
||12|| ). We are going to construct an operator function on this domain, i. e. a mapping 
F — ► <p(F) which satisfies (p(U FU^ 1 ) = Uip(F)lJ- 1 for an arbitrary isometry U. We are 
going to prove that this function has the following properties: 

<p{-F*) = -<p(F), (4.3) 
V (F + G) = <p(F) + <p(Gt), (4.4) 

for commuting F and G. This function will be built from the two gauge- invariant, self- 
adjoint and commuting (Casimir) operators K2 and K3 of F: 

K 2 = F A B F B A (4.5) 
K 3 = \ {F A B F B C F C A + F A B F C A F B C ) . (4.6) 

The Hilbert space Ti splits into the direct sum of subspaces 7i a on which F acts 
irreducibly. Each of these subspaces is a common eigenspace of K2 and K3. Denoting the 
highest weight characterizing a given irreducible representation by (m,n), with m and n 
being nonnegative integers, the eigenvalues k 2 and k 3 of K 2 and K 3 are given by: 

2 

k 2 = - (m 2 + ran + n 2 + 3m + 3n) , (4-7) 
3 

k 3 = -(m - n)(3 + 2m + n)(3 + m + 2n) . (4.8) 
9 

It is easy to check that the above formulae may be uniquely solved with respect to m and 
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n yielding: 



m = M(k 2} k 3 ) := v/|(fc 2 + 2) ( cos ( - arccos 3 :_ 

( 1 v^s \\ /x 

2 cos - arccos — = — 1 , (4.9) 

n = N(k 2 ,k 3 ) := — \/'~(k 2 + 2)[ 2 cos I - arccos ■ 



l VQh 





3 v/(A; 2 + 2)3 

cos I - arccos — 3 _ i J _ ]_ _ (4-10) 



3 v^+2; 
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Using these functions we may define a function with values in Z3, which may be identified 
with the center C of the gauge group SU(3): 

f(k 2 ,k 3 ) := (M(k 2 ,k 3 )-N(k 2 ,k 3 )) mod 3. (4.11) 

For our purposes it is convenient to use the parametrization Z3 = (—1,0,1). Since K 2 
and K 3 are commuting and self-adjoint, there exists an operator-valued function: 

<p(F)=f(K 2 (F),K 3 (F)). (4.12) 

This means that ^p(F) may take eigenvalues —1,0, 1 and that every irreducible subspace 
l-i a is an eigenspace of <p(F) with eigenvalue m — n mod 3. 

To prove property ( |4.3|) observe that K 2 {— F*) = K 2 (F), whereas K 3 (— F*) = —K 3 (F). 
Consequently, due to ( |4.7| ) and ( |4.8| ), we have M(k 2 , —k 3 ) = N(k 2 , k 3 ) and N(k 2 , —k 3 ) = 



M(k 2 , k 3 ) which implies ( |4.3j ). It remains to prove property ( [4.4| ). Let there be given two 
commuting representations, F and G, of su(3) in 7i. Denote the irreducible components 
of F by {Ti-a} an d of G by \7i G }. The irreducible spaces may be chosen in such a way 
that Ti decomposes as follows: 

W = 0W£n?#. (4.13) 

a,/3 

Take 7^ x G H F fl TL G and consider the space H F+G C 7i generated by vectors 
{F(g)G(g)x}, g G SU(3) , where F and G are the corresponding ("integrated") rep- 
resentations of SU(3). By construction, 7~C F carries an irreducible representation of 
FG . There exists a canonical embedding 

H F+G 3y^T(y)eH F ®H G , (4.14) 

given by 

T (F(g)G(g)x) := F(^)x ® , (4.15) 
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intertwining the representation FG with F ®G. This means that FG, acting on 7i^ +G , is 
equivalent to one of the irreducible components of F<g)G. Passing again to representations 
of the Lie algebra su(3), we conclude that F + G, acting on H^ +G , is equivalent to one of 
the irreducible components of F © 1 + 1 © G, acting on H^. ® TCg . Now, property (|4.4|) 
follows from the following 

Lemma 1 Let there be given two irreducible representations (m, n) and (m',n') ofsu(3), 
together with the decomposition of their tensor product into irreducible components, 

(to, n) © (to', n') = (toi, ni) © . . . © (to p , n p ) . (4.16) 

T/ien we /iat>e 

(m — n) mod 3 + (m' — n') mod 3 = (m^ — rij) mod 3 , (4-17) 

for every i — 1, . . . ,p . 

Proof: We use the classification of irreducible representations in terms of Young- 
tableaux. For su(3), the following two equations hold: 

1. The number of boxes constituting the Young-tableau of (m^, n*) equals the sum of 
boxes of the tableaux corresponding to (m,n) and (m',n ! ) minus 3p, where p is a 
nonnegative integer. 

2. The number of boxes of an arbitrary irreducible representation of su(3) is equal to 
m + 2n = m — n + 3n . 

Taking the first equation modulo three and using the second equation yields the thesis. 

Applying the operator function (p to the local Gauss law ( |3.33| ) and using additivity 
(|4.4j ) of if we obtain: 

J2v(E(x,y)) = <p(p(x)) . (4.18) 
y 

This is a gauge invariant equation for operators with eigenvalues in Z 3 , valid at every 
lattice site x. The quantity on the right hand side is the (gauge invariant) local colour 
charge density carried by the quark field. 

5 The Global Charge and the Flux Law 

Using the commutation relation between E and U, transformation law ( |2.27| ) for E(x,y) 



may be rewritten in three equivalent ways: 

E A B (y,x) = U A D (y,x)E D c (x,y)U c B (x,y) + ^6i (5.1) 

= -U c B (x, y)E D c (x, y)U A D (y, x) - p A (5.2) 

= -E D c (x,y)U A D (y,x)U c B (x,y) . (5.3) 
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These equations imply that 



K 2 (E(x,y)) = K 2 (E(y,x)) , K 3 (E(x,y)) = -K 3 (E(y,x)) . (5.4) 

Hence, we have: 

<p(E(x,y)) + <p(E(y,x)) = 0, (5.5) 

for every lattice bond (x,y). 

Now we take the sum of equations ( 4.18| ) over all lattice sites x G A. Due to the above 
identity, all terms on the left hand side cancel, except for contributions coming from the 
boundary. This way we obtain the total flux through the boundary dA of A: 

$ aA := ^<p(E(x,oo)) , (5.6) 

xedA 

where by E(x, oo) we denote the colour electric charge along the external link, connecting 
the point x on the boundary of A with the "rest of the world" . On the right hand side 
we get the (gauge invariant) global colour charge, carried by the matter field 

Q a = 5>(pO«0)- (5-7) 

xGA 

Both quantities appearing in the global Gauss law 

= Qa , (5.8) 

take values in the center Z 3 of SU(3). The "sum modulo three" is the composition law 
in Z3. 

In the above discussion we have admitted non-zero values of E(x, 00) at boundary 
points x G dA . In the remainder of this section we make some remarks on the nature of 
these objects. Our discussion will be rather heuristical, some points will be made precise 
in a subsequent paper. 

1. One may treat A as a piece of a bigger lattice A. Then the boundary flux operators 
E(x, 00) belong to J-'(A) and commute with J-'(A), (and also with O(A) ). They are 
external from the point of view of J 7 (A) and measure the "violation of the local 
Gauss law" on the boundary dA : 

E{x,oo) := p(x)-J2 E (y^) ■ ( 5 - 9 ) 

y 

Non- vanishing of this element is equivalent to gauge dependence of quantum states 
under the action of boundary gauges g(x) G SU(3) , x G dA. Let us discuss this 
point in more detail. Every irreducible representation of SU(3) is equivalent to 
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some tensor representation. More precisely, denote by T™(C 3 ) the space of m- 
contravariant, n-covariant tensors over C 3 , endowed with the natural scalar product 
induced by the scalar product on C 3 . Let T™(C 3 ) C T™(C 3 ) be the subspace of 
irreducible, i. e. completely symmetric and traceless tensors. These tensors form a 
Hilbert space T(C 3 ) defined as the direct sum 

T(C 3 ) :=0T™(C 3 ). (5.10) 

m,n 

Under gauge transformations at x G <9A, physical states of QCD on A behave like 
elements of T(C 3 ), whereas the subspaces T™(C 3 ) correspond to eigenspaces of the 
invariant operators N(E(x, oo)) and M(E(x, oo)), constructed from external fluxes 
( p.9| ). If one wants to include all these gauge invariant operators into an axiomatic 
formulation, as given in Section |2], one has different options. The remarks at the 
beginning of this point suggest to postulate that these operators commute with 
all elements of the observable algebra 0(A). This corresponds to treating external 
fluxes as purely classical objects, describing extra superselection rules (cf. ||, ||, 

2. Our results concerning the charge superselection structure of QED on a finite lattice 
|| suggest, however, a second option. In |$| all irreducible representations of the 
observable algebra were classified in terms of the global electric charge Q\ contained 
in A. Representations differing only by the local electric flux distribution over the 
boundary dA, but having the same value of the global flux $qa were proved to 
be equivalent. The redistribution of fluxes is obtained by the action of certain 
unitary operators, see |5J, acting on the quantum state under consideration. Such a 
redistribution operator has the following (heuristic) counterpart in continuum QED: 

U{n) := exp ^— J n(x) ■ A(x) d 3 x^j , (5.11) 

where n = (n k ) is a divergence-free (i. e. fulfilling dkn k = 0) vector-density on 
E C M 3 . Formally, we have: 

E k {x) := {U*{n)E{x)U{n)) k = E k {x) + n k {x) . (5.12) 

It is obvious that replacing the field E by E and leaving all other observables un- 
changed gives an equivalent representation of the observable algebra. Nevertheless, 
the flux field on the boundary <9E of the domain E is changed by n _L (x), where "_L" 
denotes the component orthogonal to <9E. 

In a subsequent paper we are going to present a similar result for lattice QCD. 
We shall prove that all irreducible representations of the observable algebra 0(A) 
of QCD on a finite lattice are classified by the value of the global colour charge 
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Q\ , yielding three different superselection sectors labeled by elements of Z 3 . How- 
ever, the local distribution of the (gauge invariant) gluon and anti-gluon fluxes 
M(E(x, oo)) and N(E(x, oo)) over the boundary dA may be arbitrarily changed 
within one sector. The redistribution of fluxes is obtained by the following pro- 
cedure. Take an arbitrary pair of points £, rj 6 dA at the boundary and a path 
(collection of lattice links) 7 = {(£,a;i), (xi,x 2 ), . . . , (xk,r/)}, connecting them. De- 
fine the following operator- valued matrix U(j) = (U A b('~y)), where 



U A b{i) 



U A c^x 1 )U c ^ c ,{x 1 ,x 2 ) . . .U c « B (x k , V ) . (5.13) 



The action of U{^) on a quantum state ip is, by definition, a collection {U A b{i)^) 
with an extra contravariant index A at £ and an extra covariant index B at r\. In 
general, the new state does not belong to any irreducible representation of SU (3) 
at £ and rj, even if ip did. This means that U^ip is not an eigenstate of operators 
M(E(£, 00)), N(E(£, 00)), M(E(r), 00)) and N(E(r], 00)), even if ip was. Decom- 
posing it into irreducible representations, we observe, however, that the value of 
<p(E(£, 00)) has been changed by plus one and the value of ip(E(r], 00)) has been 
changed by minus one by this procedure. This suggests, that these objects should 
be rather treated as quantum and not as classical quantities. Only their sum, the 
global flux c&dA , is a classical object proportional to the identity on every superse- 
lection sector. This point of view was strongly advocated by Staruszkiewicz already 
a decade ago (see fll3"lQ . We also refer to Giulini [ftlH , who discussed decoherence 



phenomena in QED in terms of quantum fluxes at infinity. 

Finally, we stress that it does not make sense to attribute any physical meaning to 
the external gluon or anti-gluon fluxes M and N themselves. It is only the quantity 
(M — N) "modulo three" which makes sense. Heuristically, this again can be made 
transparent for the continuum theory: Indeed, if we want to assign the value of the 
flux through a piece S C <9E of the boundary of a domain E, we must be able to 
integrate E over S. Suppose, for that purpose, that S has been divided into small 
portions, 5* = |J Q S a . We have 

E = J2 I E. (5.14) 

The functions M, N or even (M — N) are not additive and, therefore, cannot be 
applied to the left hand side in a way which is compatible with the Riemann sums 
arising in the integration. On the other hand, the function ip is additive. This 
enables us, in principle, to define the local Z,3-valued flux through S as the sum 
of fluxes corresponding to its small portions S a . Hence, if (one day in the future) 
continuum QCD will be constructed as an appropriate limit of lattice theories, 
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also these surface fluxes will be defined as limits of appropriate Riemann sums, 
corresponding to these lattice approximations. We conclude that the only additive, 
gauge invariant aspect of both the colour charge density p and the surface fluxes E 
is carried by the Z 3 -valued quantities ip(p) and (p(E). 
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